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ABSTRACT 



Aims. We investigate a method to increase the efficiency of a targeted exoplanet search with the transit technique by preselecting 
a subset of candidates from large catalogs of stars. Assuming spin-orbit alignment, this can be done by considering stars that have 
higher probability to be oriented nearly equator-on (inclination close to 90°). 

Methods. We use activity-rotation velocity relations for low-mass stars with a convective envelope to study the dependence of the 
position in the activity-v sin ; diagram on the stellar axis inclination. We compose a catalog of G-, K-, M-type main sequence simu- 
lated stars using isochrones, an isotropic inclination distribution and empirical relations to obtain their rotation periods and activity 
indexes. Then the activity - v sin / diagram is filled and statistics are applied to trace the areas containing the higher ratio of stars with 
inclinations above 80°. A similar statistics is applied to stars from real catalogs with log(i?J^^) and v sin i data to find their probability 
of being equator-on. 

Results. We present the method used to generate the simulated star catalog and the subsequent statistics to find the highly inclined 
stars from real catalogs using the activity-v sin ; diagram. Several catalogs from the literature are analysed and a subsample of stars 
with the highest probability of being equator-on is presented. 

Conclusions. Assuming spin-orbit alignment, the efficiency of an exoplanet transit search in the resulting subsample of probably 
highly inclined stars is estimated to be two to three times higher than with a global search with no preselection. 

Key words, stars: activity - stars: rotation - planets: detection 



1. Introduction 

In the past years, the great development in the exoplanet de- 
tection techniques and the new space missions and ground- 
based instrumentation is yielding a high rate of new discoveries. 
Transiting planets represent a treasure in exoplanet research as 
they give us the possibility to determine both their size and mass 
(the latter with a radial velocity follow-up), a nd also to study 
in de tail the properties of their atmospheres jLopez-Morale si 
I2OIO1) . The stellar light filtered through the planet's atmo- 
spher e during a transit allows to obtain the trans mission spec- 
trum dCharbonneau et alj|2002l: iTinetti et alj|2007h . whereas the 
planet's day side emission spectrum can be obtained during sec- 
ondary eclipses, yielding measurements of the composition and 
ther mal structure of the planet's atmosphere arnong o ther proper- 
ties (iBurrows et al ■l2005Ll2006LlGrillmair et alJ2008h . Exoplanet 
searches using the transit technique are nowadays providing a 
large number of new findings. Since the first observation of a 
transit for t he pla net HD 209458 b dCharbonneau et alJl200(i 
iHenrv et alJ llOOOh . 133 transiting exoplanet^ have been de- 
tected and confirmed to date, more than 1200 ca ndidates from 
the K epler mission are awaiting confirmation ( Borucki et alJ 
l2011i.b) and more discoveries are expected from several ongo- 
ing surveys. Detecting an Earth-like planet transiting a Sun-like 
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bright star is one of the main objectives of exoplanet research 
today. 

Most exoplanet transit detection programs that are currently 
underway are focused on large catalogs of stars with no pre- 
selection, basically performing photometry of every possible tar- 
get up to a certain limiting magnitude. This necessarily makes 
such surveys quite inefficient, since large amounts of data are 
processed for a relatively low transiting planet yield. However, 
some stellar properties could be used to select which stars may 
stand the best chances for finding transiting planets. This can be 
especially important for some space missions, such as CoRot and 
Kepler (and future missions like PLATO and TESS), which must 
pre-select the targets and only downlink the data from the spe- 
cific pixels containing the stars of interest. One possible way to 
perform a pre-selection is a metallicity -biased survey. In the case 
of solar-type stars it is known that planet existence is strongly 
correlated with the presence of heavy elements in the host star 
(iMarcv et al.l l2005l e.g.,). Stellar spectral type and age could 
also be considered in the case of searches of terrestrial plan- 
ets in the so called 'Habitable Zone - HZ' fealtenegger et al.l 
l20Tft) . which is defined as the region around a star where one 
would expect conditions for the existence of liquid water on a 
planet's surface. This depends on the stellar luminosity, and so 
does the period of the planets orbiting in the HZ. An example 
of a targeted sur vey with a pre-selection of low-mass stars is the 
MEarth project dCharbonneau et ani2008h . While the preselec- 
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tion increases the probability of detecting transits, limiting the 
target list to given stellar properties can result in the introduc- 
tion of selection effects on the properties of the exoplanets to be 
discovered that need to be properly accounted for 

Considering a general targeted search of transiting planets, 
little effort has been made in the past to put constraints to the 
input star catalogs with the goal of increasing the tra nsit de- 
tection rate. Relevant ideas were presented by Beatty & Seageii 
(l2010l) : transit probabilities can be enhanced if we are able to 
constrain the inclination of the stellar axes, and this can signif- 
icantly lower the number of targets to be observed in a transit 
survey. A targeted transit search would imply to observe only 
a specific sample of bright stars, which would be spread over 
the entire sky. This would require several observatories at differ- 
ent latitudes and, in principle, a significant amount of telescope 
time, as one star would be observed at a time. The detection of 
transits or even the discovery of a habitable Earth-analog with 
this approach is unrealistic. However, if stellar inclination can 
be estimated for a large sample of bright stars and expecting 
planets to orbit close to the stellar equator plane, we could se- 
lect and observe only those stars that have higher probability to 
be equator-on, and therefore lower the number of targets and 
increase the transit detection probability. In an isotropic distri- 
bution of the stellar rotation axes, only about 17% of the stars 
wo uld have spin axis inclinations above 80°. 

[Beatty & Seageil (12010) calculate how constraining stellar 
inclination affects the transit probabilities and the reduction of 
the number of targets that need to be observed for a certain num- 
ber of findings. They also discuss some ways to measure stel- 
lar inclination, but this seems to be the limiting point for the 
application of this approach. Perhaps, the most plausible possi- 
bility is to obtain stellar inclinations from spectroscopic v sin / 
measurements. For this, the true rotational period of the star 
could b e determined from photometric modulat ions caused by 
spots (M essina & G uinan 2003.' Lane et al.ll2007[) or mod ulation 
of the Ca II H and K emission fluxes ( Noves et al.lll984 ). and 
the stellar radius could be obtained from stellar models. The ma- 
jor shortcoming of this method arises from the uncertainties in 
all the ingredients: v sin /, stellar rotational period and stellar ra- 
dius. Moreover, the behavior of the sine function itself becomes 
a drawback, as it is weighed towards sin/ = 1, thus changing 
very slowly near / - 90°. All this makes necessary to know the 
observed quantities to better than 1 % accuracy if one wants to se- 
lect the stars with / > 82° (Soderblom 1985b). With the currently 
available techniques, it is not feasible to measure sin / to better 
than 10%. This approach would also be very time-intensive, re- 
quiring high resolution spectroscopy and long time-series pho- 
tometry. Therefore, it is unfeasible to obtain a relatively large 
catalog of stellar inclination data. 

Another approach is discussed in this paper, where we con- 
strain the inclination through the relation between the activ- 
ity and the projected rotational velocity of the star Section 2 
presents the principles of the approach. In Sect.|3]we describe the 
simulation of large samples of stars used to better understand this 
relation and to determine the probability for an observed star to 
have its spin axis inclination above a specific angle. This allows 
us to select the best candidates for a targeted exoplanet transit 
search. Measurements of an activity indicator, such as log(/?^^), 
and of V sin / are necessary to constrain stellar inclinations by this 
method. A catalog selection is described in Sect.|4]together with 
the implementation of the selection method and the compilation 
of a subset of stars that are expected to have inclinations close 
to 90°. The discussion included in Sect. |6] presents some com- 
ments on the applicability of this preselection method and the 



complementarity of the subsequent targeted transit search with 
the currently ongoing exoplanet searching methods. 

2. Stellar inclination from activity and rotation 

A feasible approach to estimate stellar axis incUnations with the 
cuiTently available data is to exploit the existing relation be- 
tween stelto_acU\at};aiidj2iSli2S^^ 

stars dSoderblomll 1 9851: iPizzolato et al]|2003l: iKiraga & StepienI 
'2007*). In these objects, the regime of differential rotation at 
the convective envelope plays a key role at the generation of 
the magnetic fields through the dynamo effect. These magnetic 
fields are essentially responsible for all the phenomena that are 
globally known as stellar activity, and are also thought to be the 
main rotation braking mechanism due to angular momentum loss 
through interaction with the stellar wind. Essentially, the stellar 
mass (or spectral type), which is related to the depth of the con- 
vective layer, a nd the rotatio nal period determ ine the amount of 
stellar activity (iKitchatinov & Rudiget ||1999h. This can be mea- 
sured through sever al indicators (Sod erblomll985HMallik 119981: 
Hem pelmann et a l. 1 995) . Strong evidence exists for an activity- 
rotation relation_extendmgfr(m 
dwarfs (IPizzolato et al.ll2003HKiraga & Stepie32007h . 

If an activity indicator like log(iJ^^) can be compared to 
the projected rotat ional velocity for a large sample of stars 
(iJenkins et al.ll2011i) . then those stars with / 90° (sin/ « 1) 
will have the largest v sin / values for a certain activity, covering 
a specific area on the activity-v sin / diagram for each spectral 
type. Studying the relationship between chromospheric activity 
and projected rotational velocity through such a diagram, statis- 
tics can be performed to find the area containing preferentially 
stars with / ^ 90°, even without exactly knowing the coiTela- 
tion between activity, rotation and spectral type. A large sample 
of stars is needed in order to perform a previous study of the 
activity-v sin / relations for the different spectral types. However, 
some complications arise when trying to compile log(/?^^) and 
V sin / data, since the cuiTent measurements are too scarce and 
imprecise. A solution to this is the simulation of a large sample 
of stars (Sect.O. 

The basic assumption in a selection of highly inclined stars 
for a targeted transit search lies on the alignment between the 
stellar rotation and the planet's orbital spin axis. Although from 
conservation of angular momentum we would expect the planet 
to orbit close to the stellar equator plane, recent spectroscopic 
observations during exoplanet transits have reveale d significant 
spin-o rbit missalignments for 1 of 26 Hot Jupiters ( Tria ud et aTl 
1201 Oh through the Rossiter-McLaughlin effect (Rossiter' ll924t 
McLau ghlin' ' 1 924) . These planets are thought to have formed 
far out from the star and migrated inwards. In this process, 
planet-planet scattering and Kozai oscillations due to additional 
companions could significantly affect the obliquity of the or- 
bit (Wu & Murray 2003; Rasio & Ford 1996). In spite of scarce 
statistics, spin-orbit missalignments have only been observed in 
Hot Jupiters, and assuming that a planet has formed and mi- 
grated in a disc, it is expected the majority to be in aligned orbits. 
Several multiple transiting systems have recently been found by 
Kepler (Lissaueret al. 2011 b), giving more weight to the exis- 
tence of planets with spin-orbit aligm ent. 

It is also worthwhile noting that Lanzal (1201 Oh reported ev- 
idence that Hot-Jupiters orbiting close to their stars can af- 
fect their angular momentum evolution by interaction with their 
coronal fields (Lanza 2009). This would complicate our ap- 
proach, as slightly different rotation rates would be expected for 
stars with giant planets, preventing us to predict the stellar axis 
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inclination from the same activity-v sin / distribution as for gen- 
eral stars. However, this effect has only been observed for hot 
giant planets orbiting very close to early-type stars, and the pos- 
sible induced rotation rate bias is smaller than the typical v sin / 
precisions. Therefore, no actual complexities are added. 

3. Stellar sample simulation 

3. 1 . Aims & assumptions 

The simulation of a large sample of stars containing the basic 
data for the subsequent analysis is needed if we desire to study 
accurately the correlation between chromospheric activity and 
projected rotational velocity. The resulting sample should fol- 
low the available activity and v sin / data from different existing 
catalogs, including the effects of observational eiTors and cos- 
mic dispersions for the different parameters. A simulation of this 
kind can be achieved using evolutionary models and several em- 
pirical relations. 

Among the output data of the simulation will be the stellar 
projected inclination determined from an isotropic distribution 
of the rotation axes. Different inclinations cover different areas 
in the activity-v sin / diagram, and whereas the current available 
observed data are too scarce to properly study their distribution 
and trends, the simulation presented here will allow us to accu- 
rately study how are they distributed in the diagram and perform 
a selection of observed stars from real data using the method de- 
scribed in Sect. 13. 31 The possible bias resulting from our selected 
sample will be discussed in Sect.|6] The main idea of the method 
is not to determine stellar inclinations of real stars, but to esti- 
mate the probability for each star to have higher inclination than 
a certain value. This means, under the assumption of spin-axis 
alignment, defining a subset of stars for which a transit search 
would be most efficient. 

Our study is limited to G, K, and M dwarfs (0.6 <{B-V) < 
1.6). This is mainly because we require a convective envelope 
to assume the connection between the rotation rate and the level 
of stellar activity. Therefore, we are constraining B - V so that 
the spectral range covered by our selection strategy contains all 
the stars where the activity has been observed to be scaled by 
the rotation period. This includes up to late M-type stars, so the 
same analysis can be performed there (iMohantv & Basrill2003t 
IWest & Basrl 120091: llrwin et al.l l2011 h. However, our method 
finds some limitations when applied to M-type stars from the 
fact that magnetic activity is known to rapidly saturate at a given 
level as rotation rate increases. A more extensive discussion on 
this is given in Sect. |6] 

3.2. Generation of tine sample 

The masses of the simulated stars are generated so that they fol- 
low the distribution of the Present Day Mass Function for the so- 
lar neighborhood bv lMiller & Scald (I1979i) and limited to within 
0. 15 Mq and 1 .05 Mq. These limits account for the spectral range 
we are interested to cover The upper limit is set at stars of spec- 
tral type GO. For higher mass stars the depth of the convection 
is increasingly small_and henc e stellar activity may be signifi- 
cantly less intense (iGravlll982l) . On the other hand, we exclude 
stars later than ~M4 where as will be shown in Sects. 13.31 and 
m activity-rotation pattern leaves little chance to our selection 
method considering the precision of the current v sin / data. 

Our sample is restricted to main sequence stars as magnetic 
activity behavior in evolved stars is still not well understood. 
Therefore, the stellar ages are generated considering the main 



sequence lifetime, which depends on the stellar mass and was 
calculated by fitting the terminal age main sequence points from 
the evolutionary tracks of Pietrinferni et al. ( 2004), generated us- 
ing the BaSTI web toofl We obtained the following expression 
for the main sequence lifetime t,„s'. 

log T,„, = 9.92 -3.85 log M/Mo + 2.50(logM/Mo)2- 
-1.67(logM/M0)3 ^ ' 

In fact, the upper limit on age only affects G-type stars, as the 
ages of lower-mass stars are limited by the range covered by 
the evolutio nary models iised la ter to compute stellar radius and 
color index dMarigo et al.ll2008l) . As a result, the age distribution 
for most of the spectral range is flat, covering from 0.1 to 12.6 
Gyr. The oldest lo w-mass s tars are expected to be inactive and 
very slow rotators (lBarnesl l2007). and hence they would show 
very low values of v sin /. As will be seen later in Sects. [33] and 
m such V sin / values are likely to be far from the measurement 
possibilities of current spectrographs if we want to distinguish 
and select different ranges of stellar inclinations. Therefore, the 
exclusion of old K and M-type stars does not introduce any limit 
to the optimization method that we are designing. 

The {B - V) color indi ces are derived from mass and age val- 
ues using stellar models (Marigo et al. 2008). As we are simu- 
lating samples of stars in the solar neighborhood, interstellar ab- 
sorption is negligible and no reddening is considered. The simu- 
lated sample is limited to a specific range in {B - V) between 0.6 
and 1.6. This is important as one of our goals is to study the vari- 
ations in the activity-rotation behaviour with spectral type and 
how this can affect the possibility of resolving different ranges 
of inclination. The statistics described in Sect. 13.31 will help to 
set some criteria concerning the validity of our selection method 
at different (B - V) ranges. 

Once (B - V) and age is known for a given star, the rotation 
period can be obtained using known empirical relations. Around 
the age of most of open clusters, many ob s ervatio ns c onverge 
to follow a f '''^ spin-down law. In ' BarnesI (l2003h and iBarnesI 
(2007), observations from several open clusters are used to ob- 
tain the rotation rate as a function of spectral type for F, G and 
K stars and calibrating the age dependence using the Sun. The 
age-rotation relations by Engle et al. 201 1 (in preparation) were 
obtained from stars in the spectral range covered by the simu- 
lation and which have age determinations through different in- 
direct methods, and thus represent a more suitable approach for 
our purposes. They follow an empirical expression of the form: 

P,„,(days) = Po + a ■ [Age(Gyr)]'' (2) 

for a specific (B - V). However, as will be commented i n Sect. 
|4] very similar results are obtained if the expression of IBarnesI 
(2007) is used instead. 

Three independent empirical relations of the form of expres- 
sion |2] were obtained for G-, K- and M-type dwarfs. The coef- 
ficients of the best fit are presented in Table [T] Then, the mean 
{B - V) was calculated for each of the three samples. Knowing 
the {B - V) of our simulated stars, their period was obtained by 
interpolating the results given by the three relations. A Gaussian 
error was added to the rotational periods obtained in order to 
consider a cosmic dispersion around the age-rotation relations. 
The standard deviation for this Gaussian distribution was esti- 
mated from the standard deviation of the rotational period data 
provided by Engle et al. 201 1 (in preparation), being near 1 day 
and having a slight dependence on age. 

" http ://albione. oa-teramo .inaf . i t/ 
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Table 1. Coefficients for the age-rotation relations (expression 
|2|l used in the simulation. 



Sp. Type 


Parameter 


Value 


Std. error 


ij 


n 


T 1 

Z.iol 


T TOO 




a 


6.579 


3.609 




b 


0.763 


0.253 


K 


Po 


2.158 


2.092 




a 


10.224 


2.922 




b 


0.663 


0.131 


M 


Po 


-6.624 


13.107 




a 


25.797 


15.333 




b 


0.711 


0.268 



0.08 



0.06 



0.04 



0.02 



20 



40 60 
inclination C) 



80 




0.4 0.6 
sin(i) 

Fig. 1. a) Frequency of stars in the simulated sample depending 
on the inclination of their rotation axis towards our line of sight, 
considering an isotropic distribution resulting from the simula- 
tion, b) Frequency of sin / values in the same distribution. The 
behavior of the sine function, weighed toward sin / = 1, makes 
difficult to select highly inclined stars. The analytic functions 
that describe the distributions are also plotted in both graphs 
(dashed lines). 



The stellar radius is obtained from linear interpolation of the 
evolutionary models of Marigo et al. (2008) and then the equato- 
rial rotation velocity can be calculated from the rotational period. 
An isotropic distribution of rotation axis orientations is assumed 
to generate the inclination of each simulated star, thus giving 
about 9% of the sample with ; > 85° and 17% with / > 80° 
(see Fig. [1]). An additional Gaussian dispersion component is 
added to the resulting v sin / values to include the measurement 
uncertainty. High-resolution spectroscopy and cross-correlation 
techniques are currently yielding v sin ; measurements with er- 
rors lower than ~ 0.5 km s (Ide Medeiros & Mavoi] 119991: 
[Glebocki & Gnacihski 2003'). 

Regarding stellar activity, several empirical relations exist 
that correlate Ca ii H and K emission with rotational period 



and spectral type. Due to the lack of sufficient data for an accu- 
rate study, the description of the stellar chromospheric activity 
in terms of the stellar dynamo has been diffi cult ( Donahue et aT| 
119961: iMontesinos et al.ll200lh . iNoves et all (fT984 reported ro- 
tation periods for a sample of main sequence stars and demon- 
strated that chromospheric flux scales with the Rossby number, 
Ro = Prot/Tc, where P,ot is the rotation period and is the con- 
vective overturn time near the bottom of the convection zone, 
which is an empirical function of the spectral type. This gives 
a much better correlation than activity-period, but obtaining Ro 
is complex due to the difficulties in measuring stellar rotation 
rate and the estim ation of the turnover time based o n stellar in- 
terior models (Ki m & Demarqudl996h . lNoves et alJ(il984 ) give 
an empirical calibration for the turnover time in terms of (B - V), 
considering an intermediate value for the ratio of mixing length 
to scale heigth, a - 1.9. The fit is given by the polynomial ex- 
pression: 



log Tc = 



1.362 - 0.166X + 0.025.=c2 - 5.323x^ x>0 
1.362 -0.14JC, x<0 



(3) 



where x = l-(B-V) and t^. is the turnover time in days. This is 
used in our simulation to compute the convective turnov er times, 
and th en th e Rossby numbers for each star However, as lGilmanI 
(md) and'Noyes etaD (11984 1) show, different values of a make 
stars deviate from Eq. (O, although the best correlation is found 
for a = 1.9 considering mixing-length theory models. In our 
calculated convective turnover times, a Gaussian dispersion with 
ciogT, - 0.03 dex i s added to the log Tc values to t ake this disper- 
sion into account. Mamaiek& Hillenbrand (2008) give a rela- 
tion between the chromospheric index log{R'^^) and the Rossby 
number from a larger sample, which we use to obtain the ac- 
tivity for each of our simulated stars. The calculation of the fi- 
nal logj R'^j^) also accounts f or a G aussian dispersion as the data 
in Mam aiek & Hillenbrand! (12008 ) show. The amplitude of this 
Gaussian dispersion, criog(j?^^), is scaled with Ro being always 
near 0.04 dex. 



3.3. Statistics on the simulated samples 

All the parameters from the simulated stars are generated from 
empirical relations, so the resulting sample should properly ac- 
count for the trends that real observed stars show in the activity- 

V sin / diagram. The large eiTor bars in many of the available 

V sin / measurements (being some of them only upper limits) and 
the selection biases in the observed catalogs make it difficult 
to compare our simulated sample with observations. Therefore, 
a careful selection of observational catalogs will be necessary 
to correctly check the agreement with the simulated data. The 
Gaussian dispersions introduced to the different empirical re- 
lations at several steps of the simulation and the observational 
uncertainty added to the obtained v sin / values, as described in 
Sect. 13.1! are critical to obtain a simulated sample which prop- 
erly fits real data. 

As expected, for all the spectral ranges our simulations show 
a clear increase in the mean projected rotation velocity for higher 
chromospheric activity (see Fig. |2]). Also, stars with inclina- 
tions close to 90° are placed at the right side of the distribu- 
tion, thus defining the envelope region of interest, clearly distin- 
guishable from stars with other orientations. The contamination 
comes from the cosmic dispersions and observational uncertain- 
ties. As it is evident from Fig.|2] the distribution in activity-v sin / 
becomes more densely concentrated for later spectral types, be- 
cause inactive stars tend to be very slow rotators, thus requiring 
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v-sin(i) (km/s) 



Fig. 2. A sample of 100,000 solar-type stars (0.6 < (B-V) < 0.7) 
simulated with the procedure described in Sect. 13.11 Random er- 
rors with a Gaussian distribution (cr - 0.5 km s"' ) were added to 
the V sin / values. Stars with projected axis inclinations above 80° 
are represented in red and trace the envelope region at the right 
of the distribution. For a constant log(/?^^), stars are expected 
to have a very similar equatorial velocity in terms of the rota- 
tion and activity evolution assumptions, and so different vsin/ 
are mainly due to different axis projections. The histograms in 
the two overlayed plots show the density distribution of the sim- 
ulated samples at two cuts made at constant log(R'^^) = -4.8 
and V sin / = 7.5km/s respectively. 



very precise v sin / measurements in order to apply our selection 
method. The discontinuity in the exp ression used to compute 
log(R'jjj^) from the Rossby numbers (iMamaiek & Hillenbrand 
I2OO8 ) for the simulated stars causes a small bump in the activity- 
V sin / distribution near log(/?^^) - -4.35, which is more evident 
for later spectral types (see Figs.|2]and[3]l. Although this could be 
avoided by usin g a single expression for the calculation, like the 
one obtained bv lNoves et al.l (|1984|) . it would be at the expense 
of losing precision in the general shape of the distribution. 

With sufficient stars in the activity-v sin / diagram and the 
possibility of generating a sample for any spectral type, we can 
define an efficiency parameter for each location in the diagram. 
This is defined as the probability for a star in that location to have 
an inclination angle above a certain angle, and can be calculated 
by defining a small region in the diagram around the star and 
then computing the ratio: 



Nt 



(4) 



where Ni>a is the number of simulated stars in the region with 
axis inclinations above a given a angle and Nt is the total num- 
ber of stars in the same region. For a specific (B - y)-limited 
sample, and given fixed cosmic dispersions and Gaussian errors 
for the observables, the efficiency e is a property of each loca- 
tion in the activity-v sin / diagram, and therefore regions of inter- 
est for transit searches can be studied. To quantify the increase 



effectiveness of the present methodology, we assume a - 80° 
as most of the known transiting planets are found above this in- 
clination angle. In a global sample with an isotropic distribution 
of rotation axes, we would expect 17.3% of the stars to have 
; > 80°, so that - 0.173. We may define the normalized effi- 
ciency as: 



P = 



1 Ni 



!>80» 



£80" 



67" 67" 



(5) 



so that P = 1 for a non-selected sample and P > 1 for selected 
subsamples with increased high inclination probability. 

While the efficiency for a non-selected sample of stars would 
be e a; 0.173, this increases as we move to regions at the right 
side of the distribution in the activity-v sin ; diagram, so increas- 
ing P. For simulated samples limited in spectral type, an enve- 
lope polynomial function of the form y - + a2/x was calcu- 
lated by fitting the subset of stars with / > 85°, splitting the func- 
tion into two in order to account for the inactive and the active 
part of the distribution. Successive shifts were aplied to the enve- 
lope function while the efficiency was calculated considering all 
the stars lying in the region below it (see Fig.[3]l. One can study 
the efficiency changes and test the applicability of the selection 
method for different specral types by scanning the activity-v sin / 
diagram with the envelope function that traces equal-inclination 
regions and performing such statistics. 

Four samples of different spectral types are represented in 
Fig. [3] In all cases, highly inclined stars trace the envelope re- 
gion of the distribution, but considering the observational and 
cosmic dispersions, very inactive or active stars (at the satura- 
tion zone) for later spectral types become indistinguishable in 
the activity-v sin / diagram in terms of stellar inclination. The 
empirical errors of the activity calculation (see Sect. 13. Il l and a 
cr = 0.5 km s ' for v sin / were considered in all cases. The range 
of (B - y) is also critical, as rotation rate can change significantly 
with stellar mass in the Main Sequence. Figure [3] highlights this 
and also shows how efficiency decays as the width of the enve- 
lope region increases towards the upper-left part of the distribu- 
tion and more stars are included in the statistics. 

The envelope function previously obtained is first shifted a 
certain amount to increasing vsin/ until all active stars lie at 
lower V sin /. Then, it is shifted back to decreasing v sin / at very 
small steps, considering a slope so that it scans regions of equal 
inclination. The two lines define the boundaries of the region 
of interest where the statistics is calculated. As expected, the 
efficiency decreases as more stars from the part of the distribu- 
tion with lower v sin / are included in the region. For the earlier 
spectral types, a region containing the 15% of the entire sam- 
ple yields P x 2.0, which represents doubling the probability 
of / > 80° with respect to a non-preselected sample. M-dwarfs 
show a less discriminating distribution in the diagram and thus 
the efficiency only reaches ~ 1.5 (50% increase than a non- 
preselected sample) for a small region containing 1.5% of the 
total simulated sample. 

A similar approach to the previous statistics permits the cal- 
culation of the efficiency at a specific location in the diagram. 
For a real observed star that is to say the probability for this star 
to have an inclination angle above 80°, and thus evaluate its suit- 
ability as a candidate for a transit search. Instead of performing 
the calculation of the efficiency ratio (Eq. |4|i in a wide region 
as before, a small region around the point of interest can be de- 
fined. For real data, we determine this region from the observa- 
tional uncertainties in v sin / and log(P^j^). Also, observational 
errors in (B - V) for the star are taken into account for limit- 
ing the spectral range of the simulated sample used to calculate 
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Fig. 3. Four simulated samples limited in colour index as indicated. Stars with / > 80° are plotted in red. The grey envelope function, 
calculated as described in Sect. 13.3! is shifted from the right part of the distribution to the left considering a slope so that it scans 
the high inclination region (see text). At each step the efficiency (Eq. H) with a = 80°) is calculated for the subsample below the 
envelope function. The red line in the bottom diagrams shows the evolution of e as the region limited by the function is expanded to 
the left, and the black line indicates the total fraction of stars in the region (Nsrars ~ !)■ As expected, e tends to ~ 0.17 for Nstw s ~ 1- 
The dotted grey line in the upper diagrams indicates the region where e = 0.26, thus the efficiency is increased by 50% with respect 
to a non-preselected sample. 
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P. In this case, the Gaussian errors for the observables are not 
introduced when generating the simulated data points, as they 
are already taken into account when defining the box where the 
efficiency is calculated. 

4. Data selection and analysis 

The main goal of the generation of a simulated sample of stars 
is to provide the basis to perform the statistics described in Sect. 
13.31 to real data. Currently, several catalogs exist that compile 
high-resolution spectroscopy measurements of both Ca ii H and 
K fluxes and projected rotation velocities. However, special care 
has to be taken when selecting the data and cross-matching cat- 
alogs, as diff'erent authors use different spectral resolutions and 
reduction techniques, and this results in a diversity of precisions, 
detection limits and possible selection biases. 

A first test to the high-inclination selection method is to per- 
form the efficiency statistics on the stars which already have 
known transiting planets. Although sev eral of them have been 
observed to be spin-orbit missaligned dTriaud et al.ll20ldh . the 
majority of the stars with transiting planets are expected to have 
inclinations near 90°, and thus we should obtain a high efficiency 
rate when calculating the statistics with the properly generated 
star sample for each object. It is important to notice that while 
the Rossiter-McLaughlin effect provides a measurement of the 
spin-orbit angle on the plane of the sky, there is still a compo- 
nent to be determined, and this is related to the stellar inclination 
and the planet orbit angle towards the line of sight. Therefore, the 
information given by the Rossiter-McLaughlin measurements is 
complementary to the one our statistical method can provide, 
since both angles are independent. 

Also relevant are the results from lCanto Martins et a 
showing that no significant correlations exist between chromo- 
spheric activity indicator \og{R'jjf,) and planet presence or pa- 
rameters. This is important if we want to apply the described 
activity-v sin ; statistics to these stars, using the same empirical 
approach used in Sect.[3]for the observed stars. 

Lists of vsin; measurem ents for stars with planets can be 
found in Schlaufman (■2010^ and Gonzal ez et al T (l2008h . while 
[Knutson et al. (2010) compile both log(7?J^^) and vsin; data 
from several authors. The best quality measurements of the sub- 
set of G, K and M-type stars were selected and are presented in 
Table |2l When available, the error bars of the listed parameters 
were used to constrain the generated sample and the box where 
the statistics is calculated. Otherwise, a typical box size of 0. 1 
dex in log(/?^^) and 0.5 km s ' in v sin / was adopted. For each 
star, a simulated sample of lO*" stars was generated, covering a 
range around its {B - V) value (0.04 magnitudes when no errors 
are available for photometry), and the efficiency was calculated 
with the simulated stars inside the defined box in the activity- 
V sin / diagram. 

The efficiency values presented in Table |2] were calculated 
using Eq. (|4| with a - 80°, and then normalizing by the factor 
0.173 as explained in Sect. [3] The normalized efficiency (P) is 
the factor by which the probability of having / > 80° is increased 
(compared with a non-preselected sample of stars). Because of 
the uncertainties considered for the box size where the statistics 
was calculated, very few transiting host stars show a normalized 
probability P > 2. However, the mean probability of / > 80° for 
this sample is 1.41, about 40% higher than for a general sample 
(expected to give P « 1.00). From the 21 stars sample where 
the statistics was performed, 13 show P > 1.00 and 9 show 
P > 1.50. Our strategy would select these objects, and hence 



Table 2. List of G-K-M type exoplanet parent stars with avail- 
able data. Activity-v sin / distributions of simulated samples 
where used to calculate the normalized efficiency parameter P 
(see text). This gives the probability of each star to have a rota- 
tion axis inclination above 80°, divided by the same probability 
considering a sample with an isotropic distribution of rotation 
axes. 
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Kepler-6 


0.68 


-5.005 


3.0±1.0 


3.06 


1,12 


WASP-2* 


0.84 


-4.84±0.10 


1.6±0.7 


1.61 


13 


WASP-4 


0.74 


-4.50±0.06 


2.0±1.0 


0.37 


1,14 


WASP-5 


0.66 


-4.72±0.07 


3.5±1.0 


1.70 


13 


WASP- 13 


0.60 


-5.263 


2.5±2.5 


1.00 


1,15 


WASP-19 


0.70 


-4.660 


4.0±2.0 


1.56 


1,16 


XO-1 


0.69 


-4.958 


1.11+0.67 


0.80 


1,17 


XO-2 


0.82 


-4.988 


1.3±0.3 


2.78 


1,18 


HAT-P-3 


0.87 


-4.904 


0.5±0.5 


0.67 


1,19 


HAT-P-10 


1.01 


-4.823 


0.5±0.2 


0.02 


1,20 


HAT-P-U 


1.02 


-4.584 


1.5±1.5 


0.99 


21 


HAT-P-12 


1.13 
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0.5±0.4 
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1,22 


HAT-P-13 
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-5.138 
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1,23 


HAT-P-15 


0.71 


-4.977 


2.0+0.5 
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1,24 



* Spin-orbit misalignment detected through Rossiter-McLaughlin ef- 
fect. 

Re ferences: (1) Knutson et al. (2010); (2) Schlaufman (201C^ 
(3)l Gonzal ez"eral. (2008); (4) Wright et al. (2004); (5) SatoetS 
(2005); (6) Glebocki & Gnacinski (2003); (7) Narita et al. (200"^ 
(8) Sozzetti et al. (2009); (9) Bouchy et al. ( 2008); (10) Bruntt et S} 
(2010); (11) Boruckietal. (2010); (12) Dunham et al. (201cE 
(13) ITriaud et al. (2010); (14) Gillon et al. (2009); (15) Lanp 
( I2010[) ^16) Hebbetal. (2010); (17) McCuUouah et al. { 200^ 
(18) Burke et al. ( 2007); (19) Torres et al. ( 2007); (20) Bak os et al.l 
(■2009M; (21 ) , Bakos et al. t2010.); (22 ) Hartman et al., (,200^ ; (23) 
iBakos et alj ( 120091) ; (24) lKovacs et al.l ( l2oToF ! 



a number of the transiting planets in the list could have been 
found in the subsequent targeted search. 

As can be seen in Fig.|4] where two subsets of different spec- 
tral types are shown, most of the stars lie near the right side re- 
gion of the simulated samples. Four of them are even located to 
the right of the entire distribution, and thus no statistics could 
be performed around them. These cases should be investigated 
in more detail and are probably due to underestimated v sin / er- 
ror bars. Alternatively, they could be explained by anomalously 
low values of log{R'^^), corresponding to a deep minimum of 
the activity cicle, but this is deemed quite improbable. 

The next step is the application of the formalism to catalogs 
of log{R'^fr) and v sin / measurements with the aim of selecting a 
sample with higher probability of high spin axis inclinations. In 
this case, a result of P = 2.00 or even P - 1 .50 can be considered 
interesting. For example, the preselection of a sample of stars 
having Pmean - 1 -50 would represent a 50% increase on the ratio 
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Table 3. List of analysed samples of stars with log(7?^^), v sin / 
and (B-V) data. Only G-, K- and M-type stars in the main se- 
quence are included. 



5 10 
V -3111(1) (kiii/s) 



0.8 < B-V < 1.0 




v-sin(i) (km/s) 



Fig. 4. Two subsets of transiting planet host stars (red symbols): 
solar type (top diagram) and early K-type (bottom diagram). 
Individual color indices and errors are used to constrain the 
(B - V) range of the simulated sample for each star and error 
bars were used to define the box to calculate the statistics (see 
text). Note that real data points are shown for simulated samples 
of a wider spectral range for convenience. 



of highly inclined stars, a thus a 50% increase on the efficiency 
of a planet transit search assuming a spin-orbit alignment. 

As it is evident from Fig. [3] the performance of the selection 
method that we are developing is more efficient for stars between 
certain activity thresholds. It is difficult to constrain the inclina- 
tion angle for fast rotators at the saturated activity regime, as a 
wide range of v sin / values show very similar activity levels at 
this early stage of the stellar evolution. On the other hand, inac- 
tive stars for most spectral types present rotation rates that are 
too low to resolve different stellar inclinations given the disper- 
sion of the empirical relationships and the uncertainties of obser- 
vational quantities. With regards to these, the activity-v sin / rela- 
tion makes the statistics difficult or unreliable for higher (B - V) 
indices (see Sect. 13.3b . and these should be excluded from the 
selection process in order to avoid false positives. A log(/?^^) 
threshold depending on {B - V) was calculated by considering 
a vsin/ > 1 km s ' limit and applying the envelope functions 
previously used for studying the distribution (see Sect.O. These 
functions trace the right edge of the distribution for each (B - V) 
interval, and thus their cut at vsin / - 1 km s"' can be used to 
define the inactive limit. This was done for each 0. 1 mag interval 
(0.6 < [B - V) < 1 .6) and finally a simple polynomic function 
was fitted, obtaining log(7;^^)iim = -3.45 - 1.36/(B - V). Stars 
below this inactive limit are not considered for the analysis and 
selection process. 



5. A catalog for transit surveys 

Ca II H and K flux measurements of 129 6 stars made a t Moun t 
Wilson Observatory were published by iDuncan et al.l (Il99lh . 



Data source 


N. of 


Number of stars 


with 




stars 


P> 1.5 


F> 2.0 


P>2.5 


Dc91, Gb03 & Kh09 


189 


37 


29 


16 


JklO & Kh09 


509 


209 


124 


66 


Gy06, Gb03 & Kh09 


128 


37 


28 


13 


Wr04, Gb03 & Kh09 


239 


64 


44 


32 


Hr96, Gb03 & Kh09 


99 


24 


15 


9 


LplO&Kh09 


57 


9 


4 


2 


Total 


1221 


380 


244 


138 



IT99lh 



References: Dc91 (iDuncan et alj 

( iGlebo cki & Gnacinskil l 2003h . Gv06 dGray et al.l 
Hr96 (Henrvetal. 1996), JklO dJenkins et al. ' 



Gb03 
I2006h . 
SoTlh. Kh09 



jKhar chenko & Roeser 200^, LplO dLopez-Santiaeo et al.ll20Toh 
Wr04 (Wright et al. 2004). 



We converted the S fluxes, together with their uncertainties, to 
the more s tandard activ i ty ind ex log(/?^^) with the method de- 
scribed by iNoves et al.1 (11984 . The catalog was cross-matched 
with the list of about 39000 v sin / measurements compiled 
by Glebocki & Gnaciriski (2003), rejecting those that are up- 
per limits or have large uncertainties. We further cross-matched 
the resulting catalog with th e photometry from the 2.5-ASCC 
(iKharchenko & Roeserl l20()9il to obtain color indices. Finally, 
Main Sequence G-, K- and M-ty pe stars were s electe d result- 
ing in 189 objects. Additionallv. iJenkins et al.l (1201 lb present 
chromospheric activity and rotational velocities for more than 
850 solar-typ e and subgiant stars. Hipparcos distances from 
Ivan LeeuwenI ('2007") wer e used to reject evolved or subgiant 
objects, as described bv I Jenkins et al. (201 1), an d colo r in- 
dices were obtained from Kharchenko & Roeserl (l2009l) re- 
sulting in 509 Main Sequence stars in the spectral range of 
the analysis. T h e cata lo gs with chro r nosph eric a ctivity data 
of 'Henr v et al.l (Il996l) . IWright et al.l (l2004 and iGrav et al.1 
(2006) were also cross-m atched with the rotational velocities 
by Glebo cki & Gna cinski (2003) and the color in dices from the 
2.5-ASCC. Finally, Lopez-S antiago et al.l (1201 (K) present spec- 
troscopic data, including rotational velocities and Ca ii H and K 
flux for 57 main sequence stars in the 0.6 < (B -V) < 1 .6 range. 

The complete list of catalogs of G-, K- and M-type (0.6 < 
(B - V) < 1 .6) Main Sequence stars compiled here is presented 
in Table[3] together with the total number of stars and the number 
of those found to have a high inclination normalized efficiency 
over 1.5, 2.0 and 2.5. All the results, including the efficiency 
calculated for each individual star from the input catalogs, are 
shown in Table 3 (available online). 

From the 1221 stars analysed, a subsample of 380 have a 
high inclination probability increased by 50% (P > 1 .5) or more. 
This subsample contains the stars where a planet transit search 
would be most efficient. The fact that different source catalogs 
result on a different ratio of stars with P > 1 .5 is caused by the 
existence of biases towards more active or inactive stars in the 
different catalogs and the fact that our selection method is more 
sensitive to the mid and active regime of the activity-v sin / dia- 
gram (see Sect. 13.31 ). Tighter sample selections would produce a 
higher rate of transit findings while observing fewer stars. This 
is a trade-off worth considering. For example, 244 (20%) of the 
input stars show P > 2.0 and 138 (11%) give P > 2.5, which 
means more than 0.43 probability of being equator-on (/ > 80°). 
Therefore, assuming spin-orbit alignment, at least 43% of the 
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Stars that harbor planets in this subset are expected to show tran- 
sits. In general, for stars with spin-orbit aligned planets with 
Rstaria < cos(80°) and selected with P > P' , we can predict 
that at least P' cos(80°) of them will show transits. 

6. Discussion 

As shown in Sect. |5] our strategy for targeted transiting planet 
searches results in a reduction of the initial stellar sample and 
an increase in the probability of finding transits. It is important 
to stress that we are not measuring stellar inclinations directly, 
but just putting constraints via the estimation of the statistical 
parameter e (or P). This gives the probability for each star to 
have a rotation axis inclination above 80°, and thus to be oriented 
nearly equator-on. The probability can be calculated for every 
star with log(/?^^) and v sin / measurements. A preselection of 
stars with a high value of P is expected to provide a considerably 
higher rate of transiting planets than a non-preselected sample. 
Obviously, some stars with transiting planets may not be selected 
during the process and some targets with / < 80° will inevitably 
be included in the selection, but the ratio of stars with / > 80° 
will be always higher (and even more so as we increase the value 
of P) in the selected sample than in the unselected one, which is 
the main aim of our approach. 

Several aspects should be taken into account to discuss the 
credibility of the resulting probabilities, and thus the validity of 
the subsequent selected sample to serve as the input catalog for 
high-efficiency transit searches. 

Firstly, the method for selecting high inclination stars is 
based on performing statistics on simulated samples, and thus 
depends on the accuracy of the empirical relations, distributions 
and dispersions used in the simulation. As we describe in detail 
in Sect. |3] these relations were obtained from observed data of 
stars in the solar neighborhood that show certain correlations and 
dispersions. Both the expressions that describe these correlations 
and the Gaussian dispersions that best reproduce the observed 
data were implemented in the simulation of the stellar sample, 
and hence the results are simulated samples of stars that cor- 
rectly reproduce the parameters observed in the solar neighbor- 
hood. On the other hand, a better description on both the chro- 
mospheric activity and the rotation rate dependence on mass and 
age, based on more accurate data, could help to better correlate 
the simulated samples with the observations in the activity-v sin / 
diagrams. 

Secondly, the precision of the measurements of Ca ii H and 
K flux and v sin / for each analysed sample of real stars is im- 
portant and should be considered, as very large error bars would 
make the statistics uncertain and useless. Some of the currently 
published measurements of log(7?^^) and v sin / are quite im- 
precise or do not have error determinations. Ca ii H and K flux 
usually presents variability for active stars, so several measure- 
ments taken at different epochs are needed to determine the av- 
erage chromospheric activity and its uncertainty. In this work 
we considered all the objects with available log(7?^^) measure- 
ments, applying an error box of a mean size of 0. 1 dex for the 
data with no published uncertainties. For further study and more 
accurate results, stars showing chromospheric activity variations 
or having a single log(/?^^) measurement may not be consid- 
ered or should be analysed separately. In the case of v sin / mea- 
surements, these require high-resolution spectroscopy and a very 
thorough analysis. Many of the published values are just upper 
limits, have large eiTor bars or no associated uncertainties. Since 
the rotation velocities are critical at the selection process, only 
data with the best precision was used in the analysis and results. 



Finally, the introduction of some biases and selection effects 
is evident in our approach, and we should analyse whether this 
could influence the planet detection and characteristics. As de- 
scribed in detail in Sect. 13.31 the statistics provide a better dis- 
criminating power for different values of stellar inclinations at 
the top part of the activity-v sin / diagram, and this prompted 
us not to consider the most inactive stars in order to avoid a 
high level of contamination (see Sect.|4|. Therefore, our method 
is much more sensitive for active stars and the resulting selec- 
tion will be biased towards this part of the sample. This means 
that we are rejecting part of the slow rotators of the sample, and 
hence the older stars. This bias will only have some noticeable 
effect on late-type stars {{B - V) > 1.0). For example, M-type 
stars older than ~ 1 Gyr are expected to have log(/?^^) < -4.5, 
and the activity-v sin / distribution for M stars gives no chance to 
constrain inclinations in the range below this limit. 

It is important to note that the active/young range of late type 
stars is the most unexplored in exoplanet searches. Radial veloc- 
ity surveys are forced to reject active stars that tend to show high 
rotational velocities and radial velocity jitter Also, transit pho- 
tometric surveys (especially ground-based) are likely to be inef- 
ficient for this kind of stars, due to the time-varying photomet- 
ric modulations caused by starspots. On the other hand, selected 
equator-on stars resulting from our method will have estimated 
rotation periods, and hence they are suitable stars for targeted 
searches where both the signal of the spot modulation and the 
possible transits could be detected and analysed. Therefore, a 
targeted search based on selected bright active stars expected 
to have / ~ 90° would be complementary to most exoplanet 
searches cuiTently ongoing. The observables required for our 
method, mainly resulting from high resolution spectroscopy, re- 
quire specific equipments and can be time intensive when con- 
sidering large amounts of stars. On the other hand, a single mea- 
surement is needed for each star, while a photometric monitor- 
ing or radial velocity search requires long time-series for each of 
the objects. Moreover, the same data from spectroscopic surveys 
that are required by our methodology can be useful for many 
other purposes. 

Although there is a considerable amount of log(/?^^) and 
V sin / data available nowadays (see Sect. |4|, future high resolu- 
tion spectroscopic observations with better precision may help to 
obtain even more selective samples for possible targeted transit 
searches. This strategy, besides being complementary to the cur- 
rently ongoing radial velocity and transit surveys, can be more 
efficient than a global photometric search with no preselection, 
which requires multiple photometric measurements of a large 
amount of stars to result in a relatively low rate of transit detec- 
tions. In addition, targeted observations carrying out time series 
photometry of multiple objects are possible today with the in- 
creasing number of small robotic observatories. Moreover, many 
amateur astronomers are achieving high precision photometry 
and have suitable equipments to take part in a project involving 
observations of multiple stars for a transit search. The availabil- 
ity and capabilities of amateur or small telescopes may repre- 
sent the most appropiate strategy for a targeted transit search on 
bright stars. 

7. Conclusions 

The main idea of our work was to design and caiTy out a method 
to select the best stellar candidates for a transit search from con- 
straints on their rotation axis inclination. One feasible way to do 
so with the currently available data is to make a statistical esti- 
mation of the inclinations by studying the distribution of the stars 



9 



E. Herrero et al.: Optimizing exoplanet transit searches around low-mass stars with inclination constraints 



from different spectral types in the activity-v sin / diagram. The 
need to perform a simulation of stellar properties arised from 
the lack of a large database of log(/?^^) and v sin / measurements 
with sufficient quality, and allowed us to accurately study the dis- 
tribution of stars with different inclinations in the activity-v sin / 
diagram using the statistics described in Sect. l3.3l Moreover, the 
successive steps made to design the simulation chain use the set 
of empirical relations that best describe the properties of the stel- 
lar sample in the solar neighborhood. This can also be useful to 
other fields. 

Having the possibility to obtain large simulated samples of 
stars constrained in (B - V), a relatively simple statistics was de- 
signed, so that it can be performed for every object with (B - V), 
logiR'^f,) and v sin i data. As a result of calculating the normal- 
ized efficiency P (see Sects. [3]and|4|i to the stars in the currently 
available catalogs, we proved that a preselection of about 10% 
of the initial samples can be made achieving a mean efficiency 
which is 2 to 3 times better. With the assumption of the existence 
of spin-orbit aligned planets around all the stars, this means that 
an exoplanet transit search with a 3 times higher success rate 
can be designed. In our work, the application of the approach on 
more than 1200 stars with currently available data has resulted 
on a catalog containing the most suitable sample for a transit 
search. 

In the future, larger catalogs with more precise measure- 
ments of chromospheric activity and v sin / will help to derive 
more accurate relations for the simulation of the stellar samples, 
and also to obtain more reliable results for the selection of highly 
inclined stars. On the other hand, an observing strategy consider- 
ing a targeted exoplanet transit search should be designed taking 
advantadge of the incresing availability of small robotic observa- 
tories and also photometric monitoring nano-satellites that may 
be launched in the near future and that could take profit from 
the pre-selected samples of stars resulting from the presented 
method. 
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